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Abstract: In this article, we investigate the kinetic/fluid couphng on a toy 
model. We obtain it rigorously from a hydrodynamical limit. The idea is that 
at the level of the full kinetic model, the coupling is obvious. We then investigate 
the coupling obtained when passing into the limit. We show that, especially in 
presence of a shock stuck on the interface, the coupling involves a kinetic layer 
known as the Milne problem. Due to this layer, the limit process is quite delicate 
and some blow-up techniques are needed to ensure its strong convergence. 

MSC-class: 82C40, 35L65, 82D05 

Keywords: coupling, kinetic, fluid, blow-up techniques, strong traces. Burgers 
equation, hydrodynamical limit. 

1 Introduction 

Fluid models of mechanics, like the compressible Euler system, are known to be 
not valid for gas far from thermodynamical equilibrium. This is especially the 
case for rarefied gas, or gas in the context of Physics of high energy. In this case, 
it is necessary to consider the kinetic theory which involves a new 3D variable 
called variable of velocity, e.g., the Boltzmann equation. In a numerical point 
of view, this makes the system far more costly to compute. Different strategies 
have been developed to compute with kinetic theory only in a small region where 
the gas is known to be very far from equilibrium. Standard fluid schemes are 
then used in the rest of the space. We then have the question of the coupling 
between those two regions of the two kind of models. 

These kind of coupling conditions have been widely studied, especially for 
aerospace engineering. This is the case, in particular, for the coupling of the 
Boltzmann equation and Navier-Stokes system (see, e.g., J.-F.Bourgat, P.Le Tal- 
lec, B.Perthame and Y.Qiu [51, P.Le Tallec et F.Malhnger [21], or S.Dellacherie 
[12]). Most of the methods use the thermo-equilibrium kinetic function (the 
Maxellian function) to compute the in-going fluxes at the boundary. 

On the other hand, in the case of the evaporation/condensation problem it is 
well known that a layer problem is involved. In the case of Boltzmann equation. 
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this layer problem is called Milne problem. This problem has been extensively 
studied by the team of Y.Sone and K.Aoki (see ^Q\). A.V. Bobylev, R. Grzhi- 
bovskis, and A. Heintz proposed a good approximation of it in [3. Solutions 
of this kind of equation can be classified in different categories depending on 
the boundary conditions. The complete classification was obtained by F.Golse, 
F.Coron and C.Sulem in [8J. Finally, study of such kinetic layer can be found in 
asymptotic problem linked to Semi-conductors (see P.Degond and C.Schmeiser 
[ri] . or N.Ben Abdallah, P.Degond and I.M.Gamba [2]). 

In this paper we consider a simplified model for which we will be able to 
obtain rigorously the "good" coupling by performing an hydrodynamical limit 
in one region. The model was introduced by B.Perthame and E.Tadmor [27] . 
It is a caricature of the Boltzmann equation and is given by: 

dtf + ea./ ^Mf-f, (1) 

where the "equilibrium function" A4f is defined by: 

Mf{t,X,C) = l{0<C<ti(t,a;)}-l{-ti(t,2;)<5<0} 

r+oo (2) 

with : u{t,x) = / ,f{t,x,£)d^. 



The Boltzmann collision operator is replaced by the relaxation term A4f ~ f. 
This equation is a "BGK" version of the "transport-collapse" method intro- 
duced independently by Y.Brenier [6l[2 and Y.Giga and Y.Miyakawa 14J as a 
numerical scheme for scalar conservation laws. 

Consider a function u° G L°°{R). We can define the associated equilibrium 
function f ^ M{u"{x);0 e L°°{R x K) such that: 

Mf = f 

f°{x, S,) = u"{x) for every x eR. 



In [55], P.L. Lions, B.Perthame and E.Tadmor showed that the solutions G 
L°°(M+ X M X M) of the rescaled equation U) 

dtfe+^dj,^ ^^'~ (3) 

with initial value M{vP\ ■) converge, when e goes to 0, to M{u{t,x)\^), where 
u G (K+ X K.) is solution to the Burgers equation 

dtU + d^—^Q, 
u{t = 0) = M°, 
and verifies the entropy conditions 

dt^{u) + a,i/(w) < 0, (5) 
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(6) 



for every convex function </> with associated entropy flux H verifying H'{y) = 
y(j)'(y). This equation mimics the Euler equation of the compressible gas. In 
this article, our aim is to derive a natural coupling condition between equations 
((T]) and (HI) . We choose to obtain it from the asymptotic limit when e converges 
to of the following problem 

dtfe+^dj, = a,ix)iMfe - /e), for i e e M,e e K, 

f^{t = 0) for a; e R, ^ e R, 

where the function a^{x) is 1 for a; < and 1/e for x > 0. This corresponds to 
solving the problem ([T]) in the domain a; < and the problem ^ in the domain 
a; > with a boundary coupling of the type complete transmission at a: = 0. 

In most of the numerical coupling, the in-flux entering in the kinetic part is 
defined through the equilibrium function. The coupling, in this simplified model 
would be 

Jo (1) 

f{t,0-,O=M{u{t,0+):O fore<0, 

where J-^ is the positive semi-fiux of the Engquist-Osher scheme (TH] (see also 
Degond and Jin [10] for another strategy of coupling involving a smooth tran- 
sition). This provides a plausible coupling for the the system ^ (|H). However, 
this is not the one obtained in the hydrodynamical limit. This is due to the 
possible production of a kinetic layer of size e at the interface a: = (the equiv- 
alent of the one described in the context of Boltzmann equation by the Milne 
problem) (see A.Klar jlB]). We show, in this article, that the correct coupling 
obtained at the limit of (O is the following one 



dtu = 0, for i e R+, a; £ R+, 

u{t - 0, a;) = u^{x) - /_+^ f{x, d^, for x E 



(8) 



uit,x^ 0+) 
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hj^ e/(t,0-,C)rfC, for te 



dtf + £.d^f = Mf - /, for < e R+, a: e R-, C e R, 

/(i = 0,a;,e) -/°(a:,0 for a; e R^, ^ G R, (9) 

f{t,x = 0-,C < 0) = F{t,y = 0+,C < 0) for t e R+,^ e R~, 



' S^dyF = MF - F, for t e R+, ?y e R+, ^ S R, 
Fit, y = 0+, ^ > 0) = fit, X = 0-, e > 0), for t e 

/ ^Fit,y,^)dC= ^ '"^ for^eR+,?;eR^ 



(10) 
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Note that the hmit system should verify the equaHty of the flux at the interface: 



The coupling condition in ([5]) shows that it is possible, (as suggested in ([7])), to 
define the boundary condition for the fluid equation only from the out-going flux 
of the kinetic part. It is worth noticing that this condition corresponds exactly 
to the condition of Bardos-Leroux-Nedelec W , which is the natural boundary 
condition of Dirichlet type for the Burgers equation. It is defined in the following 
way 



This result, in some sense, justifies the relevance of the Milne problem in the 
case of evaporation/condensation problems. At the numerical point of view, 
the result could seem a little bit odd. The computation of the Milne problem 
is very costly, and is certainly not needed anyway. Indeed, we show that there 
exists two types of layers depending on the boundary values on each side of the 
interface. In the case of the rarefaction layer, we show that the boundary con- 
dition coincides with the numerical condition ([7]). The layer problem provides 
a different boundary condition only in the other case of shock layer. This case 
corresponds physically to a shock coming from the fluid domain that artificially 
sticks to the interface. In some way, this means that the position of the inter- 
face has not been wisely chosen. To get an accurate description of the physics, 
the full shock should be modeled by the kinetic equation. More sophisticated 
numerical schemes use a moving interface chosen on the fly, tracking exactly the 
motion of the shocks near this interface (see Goldstein and al [29]V 

For this simplified model, first steps were done by A.Nouri, A.Omrane and 
J.P. Villa ^23], C.Bourdarias, M.Gisclon and A.Omrane [4] and M.Tidriri [321I3T]. 
However, none of those works consider the existence of the kinetic layer. Let us 
also mention that F.Golse gave a detailed study of the shock profiles of ([T]) in 
the whole hue (see [15]). 

We give in the next section the precise statements of the results. The char- 
acterization of the coupling condition in ([8|) relies on the result of strong traces 
for scalar conservation laws [35] (see Panov [25] and [2^ for new versions). 

The convergence in the kinetic layer is obtained thanks to a "blow-up" 
method first applied in the hyperbolic context in [34l [33] (see also [9] for more 
elaborate results using this technique). We present this technique in section [6?T1 
It makes use of a Liouville's type lemma presented in section [6.21 

Let us fix L > 0. From now on, we will consider initial values supported in 
^ in [—L,L]. A maximum principle implies that for every e > the solutions 
to ([To]) . (O, and ([9]) preserve this support property (see Perthame [28]). This 
condition coincides with solutions u oi ^ that are bounded by L. 

The results of this paper were announced some time ago in |36j . 
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u :— V 

if and only if 

[u^/2 — fc^/2]Sign(w — v) < for every k between u and v. 
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2 Statement of the results 

We first show the weh-posedness and stabihty of the Umit problem ([8]) ([9]) (fTO)) . 



Theorem 1 Let {f,u") G L°°(R- x [-L,L]) x L°°(R+) be such that < 
Sign(^)/° < 1 and < L. Then there exists a unique solution {f,u,F) G 
L°°(M+ X E- X [-L,L]) X L°°(R+ x IR+) x L°°(E+ x M+ x [-L,L]) to 0j 
ifP|) if verifying < Sign(^)/ < 1 and < Sign(^)F < 1. Moreover, two such 
solutions (/i, Ui, Fi), (/2, U2, F2) verify for any t G M+; 

\\fl{t) - /2(t)||Li(R-x[-L,L]) + l|wi(^) - "2(<)||li(R+) 
< ll/l ^ /2 IIli(R-x[-L,L]) + ll"l - "2lUi(R+)- 

Note that we do not claim the stability of the layer F (Indeed, as we will see 
later, it is not stable). 

Then we consider the asymptotic limit. We prove the asymptotic limit in 
two different situations. The first set of initial values considered is the following. 
Let < u+ < L and consider the functions fe° G L°°{R x [-L,L]) with < 
Sign(C)/e" < 1 verifying 



fe°{x,0>M{u+;0 for x e R- , ^ £ [-L, L], 
fe°{x,£,)> M{-u+;0 for x e R+ , £, e [-L, L]. 

The second set of initial values is the following. 

Supp/"(x, ■) c] - L, M+ - r/] for a; G M^, 
f°{x,C) ^ M{-u+;^) for x eR+ ,^ £ [-L, L]. 



(11) 



(12) 



Note that this case is more restrictive, since it is constant in the fluid region. 
We then prove the following result. 

Theorem 2 Consider initial values fiP verifying and converging weakly 
to a function /° G i°°(IR x [— L, L]) in L°°*. Let f^ be the solution to (0) with 
initial value fiP. Then the family {fe\{x<o}j — J fe '^CI{x>o}) converges weakly 
to the unique solution (/, u) to (0) (0) U0\) . with initial value (/'^|{i<o}j = 

//°rfa{.>o})- 

Consider now an initial value verifying Let f^ be the solution to 

(0) associated to this initial value. Let F^ be defined by 

FS,y,0^ fe(t,sy,0, fort GR+,t/GR+,ee [-L,L]. (13) 

Then (fe\{x<o},Ue = / /e rfCI{a;>o}: ^e) Converges to {f,u,F) solution to 0) 0) 
111U\) with initial value 

/"(a;,e) for x G M-,^ G [-L,L], 
m"(x)= J f°{x,^)d^, for xeR+. 
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Note that most of the difRcuhies of this hmit he in obtaining the boundary 
condition on the in-flow fluxes in the kinetic domain. Indeed, the boundary 
conditions on the fluid part can be obtained, in a full generality of initial con- 
ditions (see Proposition [H]) . 

Note that the initial values of type pTjl can converge only weakly. This 
means that the oscillations on the initial values cannot propagate in the non- 
linear terms of the kinetic and fluid domains. However, as we will see later, 
the kinetic layer can be destroyed in the limit. Indeed, in this case, we do not 
claim convergence of the kinetic layer. But the limit functions behave as if the 
kinetic layer was present. For the second set of initial values, things are more 
tedious. In this case, the limit solution depends on the precise structure of 
the kinetic layer. It is then crucial to show its convergence. Let us first give 
the idea why such a layer should appear. Consider the function defined on 
R+ X R+ X [-L,L] by This function verifies 

edtF, + ^dyF, = MF, - F,, for t e R+, y e M+, ^ e [-L, L], 
F,(t,2;-0,O = /s(t,a; = 0,O, for t e M+, ^ e [-i, 
Passing to the limit formally, we get 

^dyF ^ MF ~ F, for t e M+,y G K+,C e [-L,L], 
F(t,y = 0,C) = /(t,2: = 0-,C), fort eR+,ee (14) 
F{t,y^ +00, $) = M{u{t,0+);O, for t e «.+ ,£, e [-L,L]. 

The last equation is obtained formally by continuity with the fluid domain. We 
will show later that, even if the limit problem can always be formulated in terms 
of the kinetic layer, in some cases this asymptotic limit on the layer can fail. 

In order to grasp this fact, we need to get a refined study of the kinetic layer 
(or Milne problem). Consider the set 

C = I (F,5) e [0, LV2] X L\R+);0 < giO ^ ^ ^ / ^^iO ^fj • 

For every boundary data {V,g) G C, we say that F e L°°(M+ x M) with < 
Sign(^)F(y, ^) < 1 is a solution to the kinetic layer problem with given data 
{V,g) if F verifies 

' ^dyF^MF~F for y G M+, ^ G (-L, L), 
F{y^O,^>Q)^g for ^e{0,L), 

(15) 

/ ^Fiy,OdC^V for yGM+. 

We show the following proposition which classifies the solution to (fT5|) . 

Proposition 3 For any {V,g) G C, there exists a unique solution to \15\). Let 
Uoo G be such that 
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Then there are two cases. 

(a) IfV = /p ^giO '"'fi have F{y,^) — for every ^ < and 

lim F(y,-) = M{uoo,-) in L\M.). 

►oo 

M^e call this case a "relaxation layer". 

(b) IfV> j;^^giOd^, have 

lim F{y, •) = M(-Uoo, •) in L\M.). 

y-*oo 

We call this case a "shock layer". 

Let us first explain the terminology. In the first case of layer, the first moment 
/ d-S, of the given kinetic data at the interface coincides with the first mo- 
ment of the given Maxwellian at infinity. In this case the process in the layer 
corresponds only to a relaxation from the kinetic data at the interface to its 
Maxwellian. Note that in this case, the outgoing kinetic flow from the layer is 
always 0. 

In the second type of layer, the value at +00 is a Maxellian supported in 

< 0} with a flux stronger than the in-going kinetic flux at the interface. This 
corresponds to a shock front moving toward the interface. We can consider this 
kind of layer as a shock profile which "sticks" to the interface. The situations in 
the two cases of layer are very different. The first set of initial conditions (|lip in 
Theorem [5] ensures that, at the limit, the layer will be always of the relaxation 
type. This implies in particular that, in this case, the refined structure of 
the layer is not needed to get the boundary condition. This is why we can 
still consider, in this case, weak limit on the initial data. Those oscillations 
can destroy the kinetic layer, but not the boundary condition. Things are 
completely different for the second set of initial data (fT2|) . which corresponds 
to shock layer. Here, the precise structure of layer is needed to get the correct 
boundary condition. This explains why this case is more complicated. A precise 
control of the function Eg, is required. 

Note that the uniqueness in Proposition [3] implies that it does not exist any 
kinetic layer matching the Maxwellian M(uoo,'C) at y = with M(— Moo,C) at 
y ~ +00. This shows that we cannot have a classical shock profile in a half 
space. This fact was already proven by F.Golse in [TS]. This situation would 
correspond to the critic situation between (a) and (h). Thus, we cannot go in a 
continuous way from a relaxation layer to a shock layer. 

The last remark has an interesting consequence. Consider the following 
initial data (non dependent on e): 

fix,0 = M{u+;0 for x < 0,^ e [~L, L], 
= M(-u+;0 for x > 0,^ e [~L, L]. 

This initial condition verifies (dJ), so Theorem [5] shows that {fe\{x<o}j''^e = 
/ /e c^CI{x>o}) converges to the function (independent of time): 

f{t,x,^) = M{u+;^) for t e R+ , x < 0, ^ £ [-L, L], 
u{t,x) = -u+ for i e M+,x > 0,^ e [-L,L]. 
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This corresponds to a steady shock locahzed at the interface x = 0. But as 
we have seen, there is no kinetic layer matching M{u~^; •) to — u"*". Hence, 
cannot converges to . PhysicaUy, this means that at the e scahng, the shock 
is swept out from the interface. At the hmit we only get the relaxation layer. 

In view of this pathology, the stability of the limit system could seem surpris- 
ing. The stability comes from the fact that, even if the kinetic layer is not stable 
with respect to the boundary data, the boundary condition /(t, 0—, ^)l{j<o} it- 
self is stable. 

3 Preliminary results 

We begin with preliminary results on the model which will be useful later. All 
those results are fairly standard and can be found, for example, in [5S]. For 
every u G [-L, L], we define M(u, •) G L°°{[~L, L]) by: 



Lemma 4 For every regular function (j) G C°°{[~L,L]) and u G [—L,L\ we 
have: 



Let N be an integer and O be an open set of K^. Then for every function 



/ G L-^{0-L\[-L,L])) verifying < Sign(0/(-,0 < 1, we define M{f) G 
L^{0-L\[-L,L]))hy: 



We have the following property: 

Lemma 5 For every function f G L°°([— i,L]) compactly supported and such 
that < Sign(^)/(^) < 1 and every convex function (j) we have: 



M{u,^) = Sign(M)l{o<{Sigim<|u|}- 



(16) 





(17) 




Moreover, if for one strictly convex function 4> we have the equality: 




then f = Mf. 
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Proof. Notice that 



Mf~f >0 on\-L,j fiOd^i 
<0 on] J f{Od^,L[. 



Let us define 



h{o= / {Mf-m)dc. 



The function h is null at — L, it is nondecreasing on ] — L, J /(^) and non 
increasing on ] J f{^)d^,L[. From Lemma S] and (fT7)) : 



^lU'^jiOdca d^ = fjiOdc 



So /i is null at L. Hence /i(^) > on [—L. L\. Finally, integrating by parts give: 



if (j) is convex. 

Now consider one strictly convex function (p. If the inequality is an equality, 
then we have: 

rL 

r{omdi = o, 

which implies that h is null for almost every ^ G [— i,L]. Because of the 
definition of h, this implies that A4f = f for almost every ^ £ [—L, L]. □ 

Lemma 6 Let g„ e ^^([-L, L])), < Sign(^)5„ < 1 and u e L^{0) 

be such that g„ converges weakly to M{u,^) in L°°'w*. Then, the convergence 
holds strongly in Ll^^{0 x [— L, L]). 

Proof. Let be a compact set of O x [~L,L]. gn is bounded in L°° so in 
L^(ri). The convergence holds weakly in L^w. But: 

[ \M{u{x),^)\^ dxd£. <Um [ \gn\^dxd^ 
Jn Jn 

< lim / |g„p dxd(, 
Jn 

<\im f Sign(Og„(x,C)dxdC= / Sign(OM(w(x), C) 
Jn Jn 

Since |M(u,^)p = Sign(^)Af (w, ^), we conclude that the norm of gn 

converges to the norm of M{u,^). Hence gn converges strongly in L^{^1) 

and so in Ll^^{0 x [-L, L])a 
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Lemma 7 Consider two functions f,g G L^{[—L,L]) verifying the compatibility 
conditions < Sign(^)/(^) < 1 and < Sign{^)g{^) < 1. Then: 

^ l{/>g} [{MfiO - MgiO) - (/(O - gim < 0. (18) 

-L 

Moreover, if this quantity is equal to then Sign(/ — g) is constant on [—L, L]. 
Especially, in this case, Sign(A^/ — Mg) — Sign(/ — g). 

By Sign(/ — g) is constant on [—L, L], we mean that either /(^) < g{^) for all 
e e i-L, L), or fiO > for all f £ {-L, L). 

Proof. E J fd^ < Jgd^ then Mf - Mg < for every ^ e [-L,L] and the 
inequality is true. If, moreover, (jlSp is an equality, then 

l{/>9} [/(O - 3(0] = r Mf>g} [MfiO - MgiO] d^ < 0, 

J-L J-L 

and so f < g on [— L, L]. 

Assume now that J f d^ > J gd^ Then Mf - Mg > ior every ^ e [-i, L] 
and: 

l{/>g}(A^/(0 - A^5(C)) < r (MfiO - Mg{0) d^ 

-L J-L 

\fio-9mdi 

-L 

L 



< 



< 



f l[f>g}{f{0~g{0)d^+ f Mf<a}{f{0~9{0)d^ 

J-L J-L 



-L J-L 

^ i{/>g}(/(e)-.9(e))dc, 



I-L 

Hence the inequality is verified. If (|18p is an equality, then 



{.no - .9(0) d^ = / Mf>9} ifio - .9(0) 

J-L 

so / > (? on [— L,L], which ends the proof □ 



4 Study of the layer 

This section is devoted to the proof of Proposition [31 Notice that it is necessary 
for (V, g) to be in C in order to have a solution to (|15p . Indeed, if F is solution 
to p^ . then integrating ^-F(0,^) with respect to ^ we find: 

y= f\F{o,Od^^ /%5(0rfe+ r ?f^(o,Orf^> r^5(0d0 

J-L Jo J-L Jo 
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Indeed ^F(0, ^) > because of the condition < Sign,^F(?;, ^). 

As said in the introduction, in the case of a relaxation layer, the value 
F{y = 0, ^ < 0) used in the limit coupled problem is 0. Hence, in this case, the 
precise structure of the layer is not needed. In some cases the layer function Fg 
does not converge to the corresponding layer solution. But this can occur only 
in the case of the relaxation layer. And we will show that, fortunately, we still 
recover the limit problem since the coupling value f{t, 0—, ^ < 0) = is verified, 
which is the same that the value provided by the layer solution. 

In the contrary, the precise structure of the shock layer is needed to defined 
the limit problem. Fortunately, in this case the layer function F^^ converges to 
the layer solution. 

We decompose the proof into three lemmas. The first one is about the 
properties of solutions to the layer problem, the second one about the uniqueness 
and the third one about the existence of the solution. 

Lemma 8 Let {V,g) G C. We set Uoo € [0,L[ such that: 

2 

Assume that F G i°°(M+x] — L,L[) verifies < Sign(^)i^ < 1 and is solution 
to the layer problem I115\) . Then, we have the following property depending on 
the boundary data. 

(a) Relaxation layer: 

^fV^!^^g{Od^ then: 

lim F(2/,.) = M(uoo,0 inL\]-L,L[). 
Moreover F{y, ^) = for every y > and ^ g] — L, 0[. 

(b) shock layer: 
ifV> Jo^^giOd^ then: 

lim F(2/,.) = M(-Uo,,e) inLi(]-i,i[)- 

Proof. We divide the proof into several steps. 

(i) First Property of relaxation layers. We consider the case of relaxation 
layers, namely V = ^g{£,)d£,. We have 



f ^giOd^+ I mo,Od^ 

Jo J-L 



But Sign(^)F > Hence F{Q,S) = for ^ < 0. Since MF is non positive for 
^ < 0, equation (fTS]) gives for ^ < 0: 

dy{^F) = MF-F<-^-^. 
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So, for every y > 0, ^ < 0: 

But Sigii(^F) > 0, hence: 

F{y, = for y e R+, e e] - L, 0[. (19) 

(a) Preliminary result for shock layers. We consider the case of shock layer, 
namely V > 
to Lemma E]): 



namely V > Jq ^giO Multiplying equation by l{o<5} we find (thanks 



is non increasing. But 

J-L Jo 

Hence for y > 0: 

/° ^F{y, Odi>V- I ^giO d^ > 0. (20) 

J-L Jo 

(Hi) Vanishing entropy up to a subsequence. Multiplying the first equation of 
p5|) by (— and integrating it on [0,y] x {—L,L), we find 

' [\~0[MF{y,0-F{y,0]d^dy 

J-L 

< t ^'mOdi.- f eF{y,Odi 

J-L J-L 

< 2 — , 
- 3 ' 

since |F| < 1. Thanks to Lemma[Sl j {—0[-^F{y^O~F{y,£,)\d£, is nonnegative, 
so it is bounded in iy^(K+) as a function of y. Hence, there exists a sequence 
yn ^ oo such that / i-0[-^F{yn,^) - F{yn,C)] d£, converges to 0. 

(iv) Convergence of A4F(yn,-) up to a subsequence. Let us denote u{y) = 
J-L -^(^' ^) d£.. Thanks to the definition to M: 

uHyn) f\MF{y^,Od^. 



2 

Using the definition of V and the result (Hi), we find that 

-V^ £^ i[MF{y^, - F{y,„0] d^ 
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converges to when n goes to infinity. Hence, up to a subsequence, u{yn) 
converges to Uoo or to —Uoo- This means that, up to a subsequence, •) 
converges in L^{] — L, L[) to M{uoo, •) or to M(— Moo, ■), when n goes to infinity. 

(v) Convergence of F{ym ■) . Since ||i^||L°° < 1, extracting a subsequence from 
the above one if necessary, F(y„, •) converges weakly in L°°* to a function Foo 
such that < Sign(^)Foo < 1. Thanks to (Hi) and the strong convergence of 
MF{yn, ■) stated in (iv) we find at the Hmit: 



Thanks to Lemma [Sj this imphes that Foe = MFoo- Hence Foo{C) is M{uoo, ■) 
or M(— Uoo,')- Thanks to Lemma [SI this imphes that the convergence holds 
strongly in L^(] — L, L[). But thanks to (i), the limit cannot be M(— Uoo, •) (at 
least if Uoo ^ 0) in the case of relaxation layers. So in this case the limit is 
M(uoo, ■)■ In the case of shock layers, thanks to (ii), A/(uoo, •) cannot be the 
limit. Hence in this case, the limit is M(— Uoo,-)- By the uniqueness of the 
limit, the entire sequence F{yn,£,) converges strongly in L^(M) to M(uoo,-) in 
the case of relaxation layer and to M(— Uoo, •) in the case of shock layer. 

(vi) Convergence for y —> oo. Consider a monotonic function h. Multiplying 
the first equation to (jlSp by this function and integrating with respect to ^, 
wc find thanks to Lemma[5]that J £,h{£^)F{y,^) is monotonic with respect to 
y. Notice that since F is bounded by 1 and compactly supported in [—L,L], 
this function is bounded and is converging to a constant when y goes to infinity. 
Thanks to (v) the limit is J ^/i(^)M(moo, d^ in the case of relaxation layer and 
/ ^/i(^)M(— Uoo, d^ in the case of shock layer. Since every regular function is 
the difference between two nondecreasing functions, this convergence holds true 
for every regular function h. This implies that for every rj > 0, The whole family 
F{y, •) converges in VQri, L[) and in VQ — L, — to the corresponding limit 
function. But since those functions are uniformly bounded in L°°{] — L,L[), 
finally the convergence holds in !?'(] — L,L[). But the limit function is an 
equilibrium function, hence, thanks to LemmalHl the convergence holds strongly 
in L\]-L,L[).n 

Let us now consider the uniqueness of solution to the layer problem. 

Lemma 9 For every {V, g) G C there exists at most one solution F to hl5\) 
verifying F G L°°(K+x] - L, L[) with < Sign(OF < 1. 

Proof. Consider Fi , F2 two solutions to problem (fT5|) for the same condition 
values {V,g) £ C. Multiplying the difference of the first equations of (fT5|) for Fi 
and F2 by Sign(i^i — F2), and integrating in ^, we find: 



dy / S.\Fi -F2\d^= / Sign(Fi - F2)[MFi - MF2 - Fi + F2] d^. 
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Thanks to Lemma [7] this quantity is non-positive. Thanks to Lemma[51 we have: 

Jo 



v=o 



and: 



Hence: 



and: 



hm / ^\Fi- F2\d^^0. 

V-^°° J-L 

/ / Sign(Fi ~ F2)[MFi ~ MF2 - Fi + F2] d^dy = Q 

Jo J-L 





^\Fi~F2\d^ 

-L 



Thanks to Lemma[7]this imphes that 9^(Sign(Fi — F2)) = 0. Let us fix a f < 0. 
Wedenotef]^ = > F2{y,0}. Notice that since -Fa)] < 4/|^|, 

{Fi — i^2)(', is continuous and so fl^ is an open subset. Assume that it is not 
empty. Denote y > one of its elements and yo = inf{z < y\]z,y[c ^^}- For 
every z £]yo,y[ and every C G] — L,L[ we have: Fi(z,^) > ^2(2,^), hence we 
have C) > F2{z, () too. This leads to: 

^dy{Fi - F2) + (Fi - F2) > MFi - MF2 > 0, on ]yo, y[. 

Hence, we find that yo = and (Fi - F2)(0,O > (-^i - -F2)(2/, ^e^^^ > 0. This 
gives a contradiction. Hence fi^ is empty and so Fi < F2 for ^ < 0. Exchanging 
the indices gives: 

Fiiy,0^F2iy,0 2/eM+, [-i,0]. 

Hence, for every y > Q 

pL pL pL 

/ ^\Fi-F2\{y,0d£.= ^\Fi-F2\{y,0d^< ^\Fi ~ F2\{0,O d^ = 0. 

Jo J-L J-L 

This implies that (Fi — F2)(y,f) — for positive ^ too. This ends the proof.D 

Let us now show the existence of the solution to the layer problem which 
ends the proof of Proposition [S] 

Lemma 10 For every condition values {V^g) G C, there exists a solution F g 
L°°(M+x] - L,L[) with < Sign(OF < I, solution to the layer problem if73)) . 

Proof. We divide the proof into several parts. 
(i) Construction of approximated solutions. 

In order to show the existence of the solution to layer problem (fTS]) , we use the 
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method of Golse [TS]. Let {V,g) € C. We construct by induction F,, in the 
following way. We set Fi = for the case (a) and Fi ~ M{—Uoo,£,) for the case 
(b) and for n > 1: 

F„+i(2;,0 = giOe'^ + r lMF,,{z,Oe^ dz fovy>0,^e]0,Ll 

Jo s 

-MFn{z,Oe^ dz fory>0,ee]-i,0[. 

We can check that F„ verifies < Sign(f)F„ < 1: 

f ^dyFn+i - MF„ ~ F„+i y > 0, C e] - L, L[, 

(21) 

[ F„+i(y = 0,e>0)=ff Ce]0,L[. 

(ii) Convergence when n — > oo. 

Let us show by induction that for every fixed point (y,^) G M+x] — L,L[, 
the sequence {Fn{y,£,)} is non decreasing. First check that 

F2{y, - Fiiy, - 5(Oe"H{^>o} > 0. 
Assume now that i^„(?/,C) > F„_i(y,^) for every (y,^) G M+x] — L,L[. Then 

= ri[X^^„(z,0-A^i^„-i(^,C)]e^rf2>0 foreG]0,L[, y > 0, 
Jo ? 

= -[XF„(z,^)-XF„_i(z,^)]eVdz>0 for^G]-i,0[, y>0. 

By this procedure we have shown that the sequence {Fn{y,^)} is non decreas- 
ing. Moreover it is bounded by 1, hence it converges almost everywhere to a 
function F(y,^) with < Sign(^)F(y, ^) < 1. By Lebesgue's Theorem, M„(y) = 
J^]^Fn{y,^)d^ converges almost everywhere to up{y) = j'^i^ F{y, ^) d^. So 
thanks to the definition to M, MFn converges strongly to Ali^ in Lj'Qj.(R+ x] — 
L,L[). Passing to the limit in the first equation to (|2T|) shows that F is solution 
to the first equation to (fT5| . The first equation of ([2T|l shows that dy{£,Fn) is 
bounded in L°° (R+ x ] - i , L [) . Since F„ G (R+ x ] - L, L [) and L~ (] - L , X [) 
is compactly imbedded in H~^{] — L,L[), from Aubin's Theorem we find that 
the convergence holds in C°(R"'"; H^^{] — L,L[)). Hence we retrieve at the limit 
n — > oo 

FiO,^>0)^g y>0, ^G]0,L[. 

(in) Flux condition inside the layer. 

We have now to show the last condition of (|15p . Let us first consider the 
relaxation case: 
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(a) V = S,g{0 dS,. First, integrating the first equation of (fTS)) with respect to 
^ we find that 

J-L J-L 

Since Fi{y,^) = for every (y,^) G K+x] — L,Q[ and F„ is non decreasing, 
we have that F{y,^) > for (y,^) e R+x] - L,0[. But at the Hmit < 
Sign(^)F(y,0 < 1 hence F{y,^) = for £ M+x] - L,0[. So 



^9(0 = e^^(o, d^ = ei^(y, Od^ y> o. 



Consider now the shock case: 

(b) V > /q ^giO dC By induction we show that for every n > 0, F„{y,-) 
converges to M{—u^, •) in L^{] — L, L[) when y ^ oo. This is obviously true for 
n = 1 and if it is true for n then, AiFn{y, ■) converges to M{uoo, •) in L^(M) too. 
The results follows from the definition to F„+i. Integrating the first equation 
of ((2T|l with respect to ^ and using the nondecreasing property of {Fn{y, £,)} we 
find: 

dy [ ^F^iy,Odi - / iMF^^,iy,0-F^iy,0)d( 

J-L J-L 

< J ^{MFr,{y,O-Fn{y,O)d^ = 0. 
Hence using the limit at H-cx), we find for every y: 

/L pL 
^F^{y,Od^> / ^M{-u^,Od(^V. 

Passing to the limit we find: 



J JF{y,OdC = VQ>V. 



In particular Vq > J ^g{£,) d£^ hence is a layer of type (b). Thanks to Lemma 

E 

lim Fiy, •) - M{-uo, •), 

I/— ►oo 

where uq G [0, L[ verifies Vb = Uq/2. Since F > Fi, we have F{y,£^) = for 
C ^ —Woo- Hence uq < Woo- So finally Vq < V and so Vq = V. □ 

We finish this section with the following Proposition which gives a property 
of good confinement of the layers. 

Proposition 11 For any (V,g) £ C, the function F solution to verifies 

\F{y,0-Foomd£,dy<^, 
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wher 



Foo - lim F{y, ■). 

y^oc 



Proof. If V = 0, then F is identically and the result holds. If not, we have 
V = {uooY/2 > 0. From Lemma [51 

with u.^, = ±Uoo- Let u{y) — J^j^F{y,^) d^. We have 

lim u{y) — u*, 

2/— >oo 

and, thanks to Lemma [5l 

{u{y) f - = {u{y))' -2V = 2 j"^ ^ ^MFiy, " F{y, 0) d( < 0- 

Multiplying the layer equation by 2^ and integrating in ^, we find, for < y < 

00, 

2dy f env. Od(^2 r aMF{y, - Fiy, 0) = Hv))' - {u.f- 
Integrating in y gives that 

oo 

[{u{y)f - {u^f]dy < oo. 

Note that —u{y)\ = ((u*)^ — {u{y))'^) / \u^, + u{y)\, and for y large enough we 
have |u* + u{y)\ > > 0. Hence 

\u{y) ~ u^,\dy < oo. 

Using MF{y,C) = M{u{y),£), we find 



/ / \MF{y,0-F^{0\d^dy^ / \u{y)-u,\dy 

Jo J-L Jo 



< OO. 



Finally, \F{y,0 - F^i^l = Sign(C - u,)iFiy,0 - F^.{0) and 

Sign{^-u,)^dyFiy,0 
= Sign(e - u,)iMF{y,0 ~ J^oo(O) ~ Sign(C - - F^{0) 

= Sign(^ - u,){MF{y,0 ~ F^oiO) - \F{y,0 ^ ^-(Ol- 
Integrating this expression in y and ^ gives 

poo pL poo pL 

/ / \F{y,0-F^md^dy< / \MF{y,C)-F^{0\d^dy+2L^ < ^. 

Jo J-L Jo J-L 

This ends the proof. □ 
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5 Well-posedness of the limit problem 



This section is devoted to the proof of Theorem [T] which shows that the coupled 
system ([8]) (|9|) (IT0l) is well-posed and stable with respect to the initial values. 

A careful reader could be surprised that the limit problem is stable with re- 
spect to the initial conditions. Indeed, in the layer, Fi —F2 cannot be controlled 
by the initial conditions. However, notice that the coupling in the limit problem 
depends only on the values at y = of the layer. And this value is continuous 
with respect to the conditions of the problem. 

In order to obtain the existence of a solution to the limit problem, we use 
the fixed point Theorem of Schauder applied on the trace of u at the interface 
X — 0. The result relies on the stability of some quantities at x — 0. 

We first show a stability result in the kinetic domain. 

Lemma 12 Let /°, /° E L°°((-cx), 0)x] - L,L[) and gi 32 e L°°{{0,oo)x] ~ 
L, 0[). Consider fi f2 £ L°°{{0, 00) x (—00, 0)x] - L, L[) solutions for i ~ 1,2 
to 

dth + id^h = Mh - /„ for t e M+, X e M-, e e] - L, L[, 
h[t = 0, X, a = fH^, for X e M-, e e] - i, i[, 
nit,x^ 0-,e < 0) = 5,(t,e < 0) for t e M+,e e] - L,0[. 
Then, we have, for every t > 



\h{t)- h{t)\didx+ I I ^\f,{s,0^,O-f2{s,0^,O\d^ds 
Jo Jo 

< f r \fi-f2\d^dx^ f f {-0\9i{s,0-92{s.0\d^ds. 

J-ooJ-L Jo J-L 

Moreover, if /° converges strongly in L\^^ to /° and ge converges weakly to g, 
then the function /e(t,0-,C > 0) converges STRONGLY to /(t,0-,C > 0) in 
^ioc((Oj c'o)x]0, where f is the solution with initial data f'^ and boundary 
data g. 

Proof. Thanks to Lemma [71 we have 

dt f 1/1 - /2I d^ + d, f ei/i - /2I < for i e e R- . 

J-L J-L 

Integrating in (f, x) gives the first result. For the second one, first notice that, 
thanks to averaging lemmas (introduced first in [16]. See [26] for our case.), 
/ fe{t,x,£,) d£_ converges strongly in Ll^^{{0,oo) x R~). Hence Mfe converges 
strongly to Mf in Ll^^{{0,oo) x W x {—L,L)). Then the Duhamel formula 
gives for ^ > 0: 

feit, 0-, a = fei-t^. Oe-' + f e'-'Mfe{s, {s - t)^, ds, 

Jo 
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which gives the result .□ 

We now give a similar result for the coupling of the fluid part with the layer. 



Lemma 13 For any g € L°°{R+ x (0,L)) with < g{t,^) < I, and u" S 
L°°{0,oo), —L < u'^ix) < L, there exists a unique solution {u,F) to 



dtu + dx— ^ 0: for i e M+ , x e 
u{t^O,x) ^u°{x), iorxem.+ , 



u{t,x = 0+) 



BLN 



kj^ £.g{t,OdL for teM+, 



' ^dyF^MF-F, for < e IR+,y G M+,e e 
F{t, y - 0+, ^ > 0) = git, f), for teR+,^e (0, L), 



+°° v(t 0+)'^ 

^F{t,y,^)dC^ ^ '"^ forieM+,y GR+. 



Moreover, if (ui,i<i) and (u2,^2) are two such solutions associated to gi, 
(respectively g2, U2), then, for any t > 0, 

/•oo /•£ A) 

I \ui{t,x)-U2{t,x)\dx+ / \£,\\F^{s,Q+,i)~F2{s,Q+,0\d^ds 

Jo J-L 

\ul{x)~ul{x)\dx+ I I £,\g^{s,0-92{s,i)\dids. 
Jo Jo 



< 



Proof. For any g, we can construct a solution to the initial-boundary problem 
of the Burgers equation (see [24] or [17] , for instance) . This solution constructed, 
the strong trace theorem [M] gives a meaning to u{t, 0+). We can then consider 
the solution of the layer constructed in Proposition [3J 

Consider now two solution (ui,Fi) and {u2, F2) of this problem. From the 
Kruzkov's theory of the Burgers equation ([19]), and using the strong trace 
theorem, we find for any t > 

1*00 

dt / \ui{t,x) ^ U2{t,x)\dx 
Jo 

< Sign(iii(<, 0+) - U2{t, 0+))(Mi(i, 0+)2 - U2{t, 0+)2)/2. 
Using Lemma [71 we find that 



dy J jmt,y,0 - F2{t,y,0\d^ <0. 
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Integrating in y between and oo and using Proposition [3] give 
|C||Fi(t,0,0-i^2(i,0,OI^^^ 

-L 

< / ei5i(i,0-52(i,0|rfe- 1"^ m{uT{t,0.i)-M{u^{t,S,).0\di, 



L 



where 



lim F^{t,y,£,) = M{uT{t),0, 1™ F^it^y,^ = M{uT{t),0- 

y — »oo y — *oo 



Note that 

-L 

= sign«(t) - ufmuTit? - 

and 

Putting those results together gives 

dt \ui{t,x)-U2{t,x)\dx+ mFi{t,0+,O-F2{t,0+,O\d^ 

Jo J-L 

< r ^\9i{t,0~92{t,0\d^ 
Jo 

+ (Sign(ui(t,0+) - U2{t,0+))iui{t,0+f - U2{t,0+f)/2) (22) 

~ (sign(ur(o - uTm<{tf - uT{tr)/2) . (23) 

Note that the two last terms have the same absolute value. Without any loss 
of generality, we can assume ui{t,0+) < U2{t,0+). If M2(i,0+)^ < ui(t,0+)^, 
then ([22|) is negative, and the sum of (I22p and ((23|) is non positive. Otherwise, 
ui{t,0+f < U2(^,0+)^ and so W2(i,0+) > and ~ui{t,0+) < M2(i,0+). Then 
the BLN condition implies that m|/2 = / ^172 d£,, and so F2 is a relaxation layer 
(thanks to Proposition [3]). This implies that u"^ — U2, and so u'2' > uf . Hence 
([23|l is negative and the sum of ([22|) and (|23|) is non positive too. This gives the 
result. □ 

We can now prove Theorem [TJ 

Proof of Theorem d Let us first show the stability of the system. We use 
Lemma [T3l and Lemma [12] with (i=l,2) 

5^(i,e>0)-/«(t,0-,e>0), ee(0,L), i>0, 
5^(i,C<0)-F,(i,0+,e<0), ee(-i,0), t>0. 
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Adding the two estimates of those lemmas gives for any t > 

i^L noo 

\fl{t)- f2{t)\dxd^+ / \ui{t)-U2{t)\dx 



-ooJ-L Jo 
fO pL p+OD 



< 



/U nl^ r+oo 
/ \f?-f^\d^dx+ / \u^,~4\dx. 
-oo J -L JO 

This gives the stabihty and the uniqueness. 

Let us show now the existence. We fix a initial data {f'^,u'^). We denote 
H{T) the set of boundary condition on the interface for ^ < 

H = {^geL\]0,T[x]-L,Q[), -I < g{t,0 <Q)- 

Note that this set is convex. We define the function T from H to H in the 
following way. We consider the solution / on (0, T) x {—oo, 0) x (— L, L) to 
the kinetic equation with initial value /° and boundary condition g. Then we 
consider u solution to the Burgers equation in (0, T) x (0, oo) with initial value 



and boundary condition (in the sense of BLN) ^2 J ^/(t, 0— , ^) d^. Finally we 

consider F solution to the layer problem with data (/(t, 0— ,^), 0+))'^/2). 
We then set 

Tig)^^Fit,0+,^<0). 

Lemma [T2l and Lemma [T3l ensure that J- is continuous from H to H. 

Consider a sequence of function gn € H converging weakly to g £ H . From 
Lemma [T2| we get that fn{t,0—,S, > 0) converges strongly to /(t,0— ,^ > 0). 
Then, Lemma Il3l ensures that T{gn) converges strongly in H to ^(g)- Hence 
J-{H) is compact. 

So, using the classical Schauder fixed point theorem, we get the existence of 
a. g € H such that Tig) = g. The associated functions (/, F, u) is then solution 

to m, mia 



6 Asymptotic limit 

This section is devoted to the proof of Theorem[21 Let us first show the following 
proposition. It states that the solution to ([51) for x > converges to the the 
solution to ([5]) with the correct boundary condition. 

Proposition 14 Let /° G L°°(]R x [-L,L]) be such that < Sign(0/° < 1- 
Denote £ L°°(E+ x E x [— L,L]) the solution to (0) with initial value f^^ . 
Then there exists e„ ^ 0, two functions f G L°°(R x [~L,L]), f € L°°(M+ x 
R X [— L, L]) such that fe\fe converge weakly to f'^,f. Moreover the function 
u e L°°{R+ X IR+) defined by 

u{t,x) = J f{t,x,^)d£, for teR+,xeR+, 
is solution to (0). 



21 



To show this Proposition, we first give a kinetic formulation of the BLN con- 
dition (see [17] for an other kind of kinetic formulation). 

Lemma 15 Kinetic version of the BLN conditions. Consider g G L°^{Q,L) 
with < < 1 and denote v — sj 2 $5(0 o^C- have 

BLN 

u := V 

if and only if there exists h G L°°{—L,0), —1 < h(^) < and m nonnegative 
measure on (— £, L) such that 

aMiu,0-{9 + hm)=m'(.0- (24) 

Proof of the lemma. Indeed, u ^■!= v if and only if we have either u — v oi 
u < —v. In the case u = v, we get ((24| with h = and m(0 ~ J^j^ C,{M{u, Q — 
giC)) dC- The function m is equal to zero for ^ < and at ^ = i. It is increasing 
on [0,m] and decreasing on [u,L]. Hence it is nonnegative. If it < —v, we take 

HO = -l[-yiCT,o] and m(0 = /^^ C(M(w, C) - (5 + hm)dC We stiU have 
m equal to for £, < u, m{L) = and m increasing for ^ < and decreasing for 
^ > 0, so TO is still nonnegative. 

Conversely, assume that u verifies (IM)) . Noting that > 0, and integrating 
in $, we find 

mV2 > 

If u < 0, this gives u < —v. Assume that m > 0. Multiplying (|24|) by (0- and 
integrating in $, we find 

0< r e^hmd^^- f m{Od^<0. 

J-L J-L 

Hence h — 0, and integrating in ^ we have 

and so u = V. This shows that if u verifies (|24l) then u ^-.^ v. □ 



Proof of the Proposition. By weak compactness, there exists £„ — *■ such 
that /„ = fe„ converges weakly to /. Using averaging lemmas, / /„ d£_ converges 
strongly to J f d^ and so Mfn converges strongly to Mf. Passing into the limit 
in the kinetic domain gives that / verifies for x e (— cxo, 0) 

dtf + id,f = Mf-f. 

For the fluid domain, multiplying the equation by e„ shows that at the limit 

Mf{t, X, - f{t, x,0^0 for almost every t > 0, x > £_ E {-L, L). (25) 

Thanks to Lemma [51 there exists to„ > such that 

ae„{x){Mfn - fn) = d^Vfln, 
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and those measures are uniformly bounded by 

nT nOQ nL nOQ pL 

Jo J-ocJ-L J-ocJ-L 

Thus, up to a subsequence, rn„ converges to nonnegative measure m (with 
possible concentration, especially at x = 0). This gives that for x > 0: 

dtf + id.J = 9^m, t > 0, a; > 0, -L < ^ < L. 

This together with ([25| is the kinetic formulation of the Burgers equation and 
so u{t, x) — J f{t, X, £,) d^ verifies the Burgers equation for x > 0. We want now 
to recover the boundary conditions on the Burgers equation. From the limit 
equation, we get also that / e BV{'Rx',W~'^''^(R'^ x {—L,L))) and so / has a 
limit from the left and a limit from the right at each point a; £ M (in the sense 
of distribution) . Moreover 

afit, 0+, - f{t, 0-, 0) - d^m{t, 0, C), t > 0, e e (-i, L). 

The strong trace theorem ensures that f{t,0+,£,) = M(u(t, 0+), ^)) for t > 
and —L < ^ < L, and the kinetic formulation of the BLN conditions ((24|) gives 

u{t, 0+) ^■.= ^ £ ^f{t,o-,0dc 

□ 

We now show the Theorem[5]in the case of the first set of initial values. We 
first consider the special initial value, for uq > 

V'°(a;,0 = l{,<o}M(uo,0 + l{r.>o}A^(-"o,e)- 

We want to show that for this initial value, at the limit e ^ we get f{t) = 
for t > 0. For this matter, we first consider the initial value 

i:°{x,o = i'°{x-s,0- 

Consider the limit function fs. We have seen that us — J fs d^ verifies the 
Burgers equation for a; > with initial value Ug{x) = l{x<s}Uo — 1[x>s}Uq- 
Thanks to the finite speed of propagation there exists a finite time (at least 
S/uq) such that u{t,0+) > 0. Hence, the BLN conditions ensures that on this 
lapse of time us{t,0+)'^ /2 — Jq (.fsit^O—, £,) d^. But the conservation of flux 
at the interface ensures then that fs{t,0—,£, < 0) = for this lapse of time 
and ^ < 0. And so, solving the equation on the kinetic domain shows that 
the solution is stationary. The values at the interface does not change on time 
either, and so the solution of the fluid domain is constant also. Finally we find 
f{t) = "05 for all t < S/uq. iterating the argument shows that this is true for all 
time t > 0. Using Lemma [71 we can show that for any e and t > 

/ \fS)- fe{t)\didx< / / -i,°,\didx = 25u^. 

-oo^—L J — oqJ — L 
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up to subsequence, the functions converges strongly, and so 

/OO nL 
/ \f{t)~i^°\d^dx<25uo. 
-OQ J —L 

Passing into the limit when 5 ^ Q gives that f{t) — ij/^ for every t > Q. 

Then, using Lemma[71 we can show that if /° < /2 then for any e, /i^^ < /2,e 
(see [55] )■ And so, at the limit /i < f2- The hypothesis on the initial data are 
exactly /° > V'^) hence, at the limit: 

f{t, X, C) > i>°{x, for t > 0, x e R, C G L). 

This implies that /(i, 0—, ^ < 0) = for t > 0, ^ < 0. But the condition implies 
also that the kinetic layer has to be of the relaxation type. The condition is 
then the good one. 

We are now left to show the asymptotic limit for the second set of initial 
conditions. We need, in this case, to pass into the limit in the term F^. To avoid 
the possibility for the layer to be swept away, we need to show that the shock is 
well confined against the interface x — Q. This will be provided by Proposition 
[TTl The hypothesis on the initial data (|12p . and a comparison argument gives 
that, uniformly with respect to e 

\F,{t,y,0-M{-u+-0\ < \Fgiy,0-M{-u+;0\ for t e M+, y e M+, ^ G M, 

where Fg is the solution to the layer problem (fT5)) associated to {g, V) with 
g{^) = M{u+ - ij; ^) and V = {u+f/2. 

We need now to get strong compactness to pass into the limit in the layer. 
For this, we use a blow-up technique first introduced in [33] . 

6.1 The "blow up" method 

The equation on F^ does not control the oscillations in time when e ^ 0. The 
idea is to get back the balanced structure of ^ doing a zoom in time of the 
equation. We introduce a new local variable of time s G K and a new rescaled 
function defined by 

F,{t,s,y,O^F,{t + es,y,i). (26) 
For almost every fixed t > 0, the function F^ verifies 

dsF, + ^dyF, = MF, ~ F,, for s g] - t/e, +«)[, y G M+, ^ e L], 

F,{t, 3,0,0 = fe{t + es,0,^) for s e] - t/e, +^[,^ e [-L, L], 

\F,it,s,y,0-M{~u+;0\ < \Fa{y,0 ^ M{-u+;0\ 

for s e]- t/e, +oo[,y eR+,£, e[-L,L]. 

(27) 

Results obtained on the rescaled functions can be translated on the non-rescaled 
one (and vice versa) thanks to the following lemma whose proof can be found 
in [33l : 
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Lemma 16 (From local to global) Let F^^, F € L°°{W x W), p,q integers. 
Then, F^^ converges strongly to F in Ll^^{MP x M'^) if and only if for every 
Ri>0, R2>0 and R3 > 0; 

[ [ [ \F,^{x + eny,z)- F{x,z)\ dydzdx""^^ 0. {28) 

JBj,{0,Ri) JBg{()M.2) JBp{()M.3) 

Let us first apply this lemma on /e(-,0, •) with p = 1 (time variable) and q — 1 
variable). From Lemma [T^ the function f^{t,0,^ > 0) converges strongly to 
fit, 0, ^ > 0) in Lioc(R+ X (0, L)). Hence for any T > 0, i? > 

\F,{t + es, 0, - f{t, 0, 0\d^ ds] dt'~^ 0. 




Hence, there is a subsequence e„ such that for almost every fixed t > 0, the local 
function (restricted to positive ^) Fe^{t, -,0, ■) converges strongly in ^^^^(R x 

(0,L))to/(t,0+,-). 

Let us fix such a time t > 0. We pass into the limit when e„ goes to 0. For 
this we first show the following proposition. 

Proposition 17 Let s„ be a sequence converging to —00 and consider Fn G 
L°°{{sn, 00) X X {—L,L)), with < Sign(^)F„ <1 a sequence of solutions 
to 

dsF^ + ^dyF^ = MFn ~ F„ , (29) 

|F„(s,y,0-M(-^+;C)| < \Fg{y,0-M{~u+-0\ (30) 
for s e]s„,+cx}[, ?; e M+,C e [-L, L]. 

Assume, in addition that i^„(-,0, •) (restricted to positive S,) converges strongly 
in LI^^(Kx (0, L)). Then, up to a subsequence, Fn converges strongly in ^[^^^(Mx 
M+ X {—L,L)) to a solution Fi of S29\) S30\}. In addition, F„(-,0,-) converges 
strongly in Ll^^{R x (-L,L)) to F;(-,0,-). 

Proof of the proposition. The function F„ is bounded by 1, so, up to a 

subsequence, it converges weakly to a function denoted Fi. We can now use the 
averaging lemmas on the kinetic equation ([^5]) . We deduce that A4Fn converges 
to MFi in LJqj,. Hence Fi verifies the same equation ([29t . Moreover, there exists 
a subsequence such that MFn and the restriction to ^ > of the trace F„(-, 0, •) 
converges almost everywhere. For ^ > 0, integrating the equation (I29|) along 
characteristics, we find that for almost any s, y and n big enough: 

F„(s, y, C) = e-^Fn{s ~ ^,0, + ^MFn{s + ^,z, dz. 

From the convergence of Fn almost everywhere on the trace {y — 0, £, > 0) and 
the convergence of MFn almost everywhere, we deduce the convergence for 
almost every (s, y, ^) G M x R+ x (0, L) of Fn{s, y, ^). This function is bounded 
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by 1. The Lebesgue's Theorem ensures the strong convergence of (the restriction 
for ^ > of) F„ in L^^iR x K+ x (0, L)). 

For ^ < 0, we use the same strategy, using ([3D]). Thanks to Lemma [51 for 
any 77 > 0, there exists y > (as big as we wish) such that for any y > Y we 
have 







\Fg{y,0-M{-u+;0\d^<V. 



Thanks to (|30p . we deduce that for any y > Y, we have for any n, m 

j sup \F^{s, y, - F^{s\ y, 0\dC < 2r]. 

Integrating, again, the equation along the characteristics, we find for n big 
enough 

— v-y— j-y eM-y'^/^ — z-v 

Fn{s,y,0 = e—F,,{s + iY-y)/^,Y,0 + — MF^{s+ z,0 dz. 

So, for n, m big enough, we have for every s, and y <Y 



J jFnis,y,0~F^{s,y,0\d^ 

/■o fv p{^-y)li _ _ z-v 
<2r,+ / \MFn~MF^\{s + —^,z,Odzd^. 



using the convergence of A^F„ almost everywhere, and the Lebesgue's Theorem, 
we find that for almost every s e M, and every y <Y ^ the restriction for ^ < of 
Fn is Cauchy in L\^^{—L,Q). Using the Lebesgue's Theorem, and the previous 
result for ^ > 0, we get the convergence of F„ in Lj'jjj.(M x x (— L, L)). 
The result with y = i) gives also the strong convergence of the trace. □ 

For almost every fixed i > 0, we apply this proposition on F^^{t, ■,-,■). 
Passing into the limit in (|27|) . we find that, up to a subsequence, F^^{t, ■) 
converges strongly to Foo{t, •, •, •) solution to 

dsFoo + ^dyFoo = MFoo - Too, for s G K, y e R+, ^ G L], 

F^{t,s,0+,O = f{t,0+,0 for seM,^e[0,L], 

(31) 

\Foo{t,s,y,0-M{^u+;0\ < \Fg{y,0 - M{-u+;Ol 
for seR,y e M+,^ G [-L,L]. 

The following Liouville lemma (proven below) completely characterizes the limit. 

Lemma 18 There exists a unique solution Foc{t, ■, ■, ■) to iSl]) . This solution 
does not depend on s. 
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Wc then deduce that Foo is the unique solution F to the kinetic layer problem 
(fTOl) associated to {f{t,0+,^ > 0),^ = (u+)2/2). By the uniqueness of the 
limit, the whole sequence is converging in Ll^^. 

Lemma [TC] ensures the convergence of the global functions ^F^^^ from the 
result of convergence of the local functions F^^ . 



6.2 Liouville's lemma 

This subsection is dedicated to the proof of Lemma [T51 Existence of a (steady) 
solution to (|3ip is given by Proposition [3] and a comparison principle. Assume 
that there exists two such solutions Fi and F2. We first have: 

ds\Fi - F2I + (dy\Fi - F2\ = F2), for s e M, y e M+, ^ £ M, 

\Fi- F2\{s,0+,C) = for seM,CeM+, (32) 

\Fi- F2\is,y,0 <2\Fgiy,0 - Mi-u+,0\ for s £ K, y £ K+, ^ £ K, 

where 

/L pL 
V{Fi,F2) = - / Sign(Fi - F2)[MFi - MF2 - Fi + F2] > 0, 
-L J-L 

thanks to Lemma [71 Thanks to Proposition [TT| Fg — M{—u^; •) is integrable, 
hence |_Fi — •? ') is integrable for every fixed s. Integrating the first equation 
in ([5^ with respect to y and ^, we get: 

-r / \Fi~F2\{s,y,i)didy 
ds J-L Jo 

/ V{FuF2)is,y,0d^dy+ ^ - F2\{s,0+,0 d^ < 0. 

-L Jo J -00 

(33) 

We can deduce that Jp°° l-^i ~ -^"2 1 (s, y, ^) is a non increasing bounded 
function. We denote 5 its limit at —00. Consider now the functions translated 
in time F[^{s, •, •) = Fi{s — n, •, •) for i — 1,2. Thanks to Proposition [T71 up to 
a subsequence, those functions converges in Lj^^ to F^,F^ solutions to (PT|) . 
In particular, the quantity \F^ — F^\ verifies (|32]) . Thanks to the Lebesgue's 
dominated convergence theorem, we get — F^ \ (s, y, ^) d^dy = S 

for every s G M. Hence, we have 



r r V{F^,Fnd^dy = 0, 

Jo J-L 



and 

\Fr -F^\{s,0,O = 0. 

The proof now follows the lines of the proof of Lemma [SI Thanks to Lemma 
[3 we have a^(Sign(Ff° - F^')) = 0. Let us fix s e R, y G M+ and ^ < 0. We 
denote 

F,iT)^Fris + T,y + T^,0, ^ = 1,2, tGR, 
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and 

^s,y,i = {0 < T < -y/e\Fl(T) > F2(t)}. 

Note that this set is open since Fi — F2 is Lipschitz (and so continuous) . Assume 
that it is not empty. Denote r one of its elements and tq = sup{A > r \ (r, A) C 
For every A £ (t, tq) and every ( £ {--L,L) we have: F^{s + X,y + 
A^! i) > F^{s + X,y + 0, hence we have F^{s + X,y + X^ C) > -^^2^(5 + 
X,y + A^, C) too. This leads to: 

dr{Fi - F2) + (#1 - F2) > MFi - MF2 > 0, on (r, tq). 

Hence, we find that To = -y/^and (Ff°-F2°°)(s-2//^, 0, > (A- A)(T)e2'/« > 
0. This gives a contradiction. Hence fi^ is empty and so F^ < F^ for ^ < 0. 
Exchanging the indices gives: 

F^{s,y,O^F^{s,y,0 seM,2/eM+, $e(-L,0). 

Hence, for every y > 



r r \F^ -F^\{s,z,Od^dz 

J-L Jo 

I ^\F^ -Fi°\{s,y,Od^^D{s,y)<0. 
Jo 



ds 



(34) 



So for any y > 0, /j'"^ \F^ — F^ \ (s, z, ^) is a non increasing bounded 
function. We denote d{y) its limit for s goes to —00. Consider a new translation 
in time 

Fr''{s,;-)=Fris-n,;-), 

for i — 1,2. As before, up to a subsequence, those functions converge in Lj^^ to 
5^2 • -Lney verity 

I 77100,00 TTiOO-OOl/ /- ^ n\ n 

\Fi -F2 |(s,y,C < 0) = 0, 

and for any y > 0, and s G K J^j^Jq \F^"°^' — F^'°°\{s, z,^) d^dz — d{y). 
Hence, from ([34)) . for any y > 0, s G R: 

Jo 

Finally F^'°° = F^'°° and d(j/) ^ for any j/. We can deduce that |Ff° - 
F^\{s,y,^) = everywhere, and so 5 = 0. This gives Fi = F2.n 
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